Quantum Error Correction in Correlated Quantum Noise 
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We consider quantum error correction of quantum noise that is created by a local interaction of 
qubits with a common bosonic bath. The possible exchange of bath bosons between qubits gives 
rise to spatial and temporal correlations in the noise. We find that these kind of noise correlations 
have a strong negative impact on quantum error-correction. 
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The superiority of quantum computation over conven- 
tional computation relies on the fact that a quantum- 
bit (qubit) register can be in the superposition of a very 
large number of classical computational states. At the 
same time, maintaining coherence of this highly super- 
positional state poses also the main obstacle for the re- 
alization of a quantum computer. 

For a small number of qubits this difficulty can be over- 
come by simply reducing the coupling to environmental 
degrees of freedom, as has been demonstrated by several 
groups for different physical implementations. With an 
increasing number of qubits, it will, however, become ex- 
tremely difficult to reach the required coherence in that 
way [lj. It is therefore common opinion that a scalable 
implementation of a quantum computer must use some 
error correction scheme that recovers the quantum state 
after it has been distorted by external noise. 

The existence of error correcting schemes for quantum 
states, which was shown independently by Shor |2( and 
Steane 0], is a remarkable fact and has been crucial for 
the development of the field. The key ideas presented 
in their work rapidly evolved to a beautiful theory of 
quantum error correcting codes and subsequently to the 
concept of fault tolerant quantum computation 

Quantum error correcting schemes are usually designed 
for the independent error model, which by definition does 
not exhibit correlations between noise of different times 
and locations. From a physical point of view, this re- 
quirement is rather annoying, since in general qubits 
do interact with a common environment which necessar- 
ily introduces some amount of correlations in the noise. 
To be more specific, in many if not all situations the 
qubits weakly interact with a common thermal bath of 
extended bosons (photons and/or phonons). The ex- 
change of bosons between qubits will then cause spatial 
and temporal error correlations that violate the condition 
of error independence. Indeed, it has been shown [5] that 
these kind of processes can lead to drastically enlarged 
or reduced decoherence of certain states. 

To which extent do such error correlations inter- 
fere with quantum error correction? We have ana- 
lyzed this problem for optimal Calderbank-Shor-Steane 
(CSS) quantum error-correcting codes of variable length 
n (number of physical qubits) and size k (number of log- 



ical qubits). As physical noise- model, we use a reduced 
spin-boson model consisting of n spins - describing an n- 
qubit register - coupled to a common bosonic bath 0, [|J . 
The amount of noise correlations can be controlled by the 
inter-spin distance r. Within this framework, we study 
how code states transform during spin-boson interaction 
and a subsequent error-correcting operation. The dis- 
tance between the resulting code state and the initial one 
- in the sequel denoted as residual error A [cf. Eq. (|ll|l ] - 
serves as a measure for the error-correcting performance 
of the code. 

Our main finding is that quantum error correction 
with CSS codes is substantially hampered by the kind of 
noise-correlations captured in our model. This becomes 
evident by the fact that for any fixed information rate 
k/n > the residual error A approaches a finite con- 
stant in the limit n — > oo, unless the spin-boson coupling 
strictly vanishes or r is infinite (cf. Fig. ^) . Using a sim- 
ple scaling argument we conclude that for a wide range 
of model parameters CSS codes cannot provide the accu- 
racy needed for large scale quantum computations. 

In contrast to related studies Jjj, here the spin-boson 
coupling is treated in a non-perturbative manner, which 
we find to be indispensable in the large n limit. Substan- 
tial progress towards error correction beyond the inde- 
endent error model has been made in very recent work 
0. We will briefly comment our results in light of 
this new work at the end of this Letter. An extended 
discussion will follow in a future publication |lf| . 

We begin with the physical model for the n-qubit reg- 
ister. It is defined by the Hamiltonian 
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of n spin-i particles (qubits) at positions r/ interacting 
with three-dimensional bosonic modes with creation (an- 
nihilation) operators b\. (6k) and energies c^k = |k| [llj . 
a z .i denotes the Pauli spin o~ z operator acting on the ith 
spin, and e is the Zeeman energy. The spin-boson cou- 
plings <7|k| may be characterized as usual by a spectral 
function J(uj) := £ k S{uJ k - u)\g w | 2 = Au'e'^", with 
a cut-off frequency f2, a constant A of appropriate di- 
mension, and a positive parameter s [l^ . Note that the 
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spin-boson couplings do not lead to energy dissipation 
but to dephasing of the spin system. 

The spin-boson model given by Hamiltonian H is quite 
suitable for our purposes because it shows full decoher- 
ence and is still analytically solvable 0, . We assume 
that initially the system is in a product state p(0) ® p% 
of a spin state p(0) and a thermal state Pfipf the bosonic 
bath at temperature T. Following Refs. pi la] we deter- 
mine the reduced density matrix p(r ) of the spin system 
at some subsequent time r to be [10( 



^( T )=E c7 n e" c - p(0) P M , 



(1) 



where P v for n 6 Z 2 = {0, 1}™ = {|, |}™ is the projector 
on the state |?7) = 1 770 ) • ■ ■ \r)n-i)i and C w are real 
coefficients given by 

C VH = E;™^ ~ W)(l7m - Mm) r |r,-r m | • 

The distance-dependent decoherence-parameter |lfjf 
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is a positive, monotonously decaying function of r if s < 
2. Up to an singularity at r = r, the same holds for 
2 < s < 3 if temperature is high, T > r _1 ,r -1 [ic| . 
We confine our considerations to s < 3 and note that 
in quantum optical settings typically s = 1, the Ohmic 
case, where IV (r) is linear in r for large t ^> 1/T [l2T |. 

To make the model manageable we simplify the coef- 
ficients C w by setting all distances |r/ — r m |;^ m to the 
maximum distance r. Keeping in mind that this simpli- 
fication lowers the effects of correlations, we obtain 



c m = h©Ai|(r 



iY) + (M 



\p\) 2 T r 
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where \rj\ denotes the Hamming weig ht H3 of ?y, and © 
means a bitwise addition modulo 2 in ZJ. 

Although representation Eq. looks simple, it turns 
out to be rather cumbersome for our further calculations. 
A much better one can be given in terms of operators Z v 
that are defined for v = (V ■ • ■ f«-i) G Z 2 as product of 
exactly those o~ Zt i where v\ = \. By the structure of the 
interaction Hamiltonian it is evident that then 



p(t) = E 



p(0) z„ 



(4) 



with suitable coefficients q„„/ (r) . To determine them we 
let p(0) = be the projector on the totally x po- 

larized spin state and compare the expressions obtained 
in both representations. In this way we see that <x vv i and 
e~ Cr *v are related by a Fourier transformation in Z 2 x Zj , 
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(i/ • 7] denotes the standard inner product in Z 2 )- Later 
on we will only need the diagonal coefficients ct vv , which 
solely depend on the Hamming weight \v\. For C w ac- 
cording to Eq. © the sum © can be easily carried out 
if we eliminate the squared term in the exponent by the 
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some algebra we arrive at 
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where we introduced the x-dependent "probability" 

Px = (f - e - r ° +rv cos2a;)/2 € [0, f] . (7) 

Thus, we have a convenient description of the decoher- 
ence process at hand and can now turn to CSS codes. 

The key idea of quantum error-correction is to encode 
the information of k logical qubits in an appropriate sub- 
space C of the Hilbert space TC n associated with n > k 
physical qubits. C C TL n is called a quantum code of 
length n and size k. Error operations that unitarily map 
C on different cosets of C can then be detected and cor- 
rected. A CSS quantum code is constructed on two linear 
classical codes [l4| (i.e. Z 2 -linear spaces) C'2 C C\ C Z2 
by 



C = span{ I Q)} Q 



where C,- is the orthogonal space of Cj with respect to 
the inner product in Z 2 , and the state vectors \Q) are 



\Q) = 




z q \y) , q^Q 



C encodes k = log 2 \C% /C'jH = dimCi — dimC2 logical 
qubits in n physical qubits 0, . The error correcting 
capability of C is determined by the minimum weig htshj 
d% and of C% resp. . An error correction scheme 
using C can correct up to t = [^^] universal qubit er- 
rors, where d — min{(ii, dj^}. This characterizes C as an 
[n, k, d] code. With P being the projector on C, the error 
correcting operation 1Z associated with C reads 



f,f*eZ5,M,|/t|<t 



PXpZvpZvX^P (8) 



(A^ is analogously defined as Z v with cr^; instead of <J z ,i)- 
CSS codes have been used to demonstrate the existence 

of efficient ( "good" ) quantum error correcting codes 0] , 

as specified by a theorem of Calderbank and Shor: 

Theorem (Calderbank and Shor): For sufficiently large 

n there exists always an [n, k, d] CSS code satisfying 



k/n > 1 - 2H 2 (d/n) = R css (d/n) 



(9) 



(H2(x) — —x\og 2 x — (1 — a;)log 2 (l — x) is the binary 
entropy function.) Provided that the error-correcting op- 
eration 1Z for such codes is error-free, this theorem can 
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FIG. 1: Residual error [A]™ of [n, k, d = It + 1] CSS codes as 
function of t for r = 0.01, k/n = 0.062, and t/n = 0.05. The 
values of F r are .01, .005, .0025, .00125 (squares, top to down), 
and F r = (dots). The dashed lines are the asymptotic values 
of the residual errors. 



be rephrased in a more pragmatic fashion: An n-qubit 
register that is perturbed in at most t qubits can be used 
to perfectly restore k = [nR css (^p)] logical qubits that 
have been encoded in an appropriate [n,k,2t + 1] CSS 
code. 

Of course, a real and noisy physical n-qubit register is 
unlikely to guarantee the prerequisite of the theorem. It 
therefore does not necessarily provide a practical solution 
of the decoherence problem. One has to demonstrate that 
also under more realistic conditions the performance of 
the error-correcting code is still sufficient. This has been 
shown in refs. 15] for the independent error model. We 
will now investigate how CSS error correction performs 
when the errors are correlated. 

To this end we consider an arbitrary state = 
St; V'»)I 7 ?)l °f an abstract fc-qubit register with 2 fc or- 
thonormal logical basis states \t])l, "q G Z§. The logical 
state \^}l is encoded in an [n, k, d] CSS code C as a vec- 
tor = J2 V tPvlQv) C according to some bijective 
mapping i— » \Qr/) ■ The encoded state = I^X^I 
is first subjected to the noise operation defined by Eq. 
(0J and then corrected by 1Z. This results in a final state 
p'q = TZ{p^,(r)) which we compare with the original state 
p^, by the fidelity F{p^, 1 p'^,) 1 which here is 

\n\<t vv' 

Since the square root of 1 — F defines a proper distance 
measure for density matrices |lfij . we call 



A* = l-Fipv,^) 



(11) 



the residual error after CSS error correction. 

Here we are not interested in specific codes but in the 
general properties of CSS codes. We therefore continue 
by taking the average over practically all good CSS codes 
of a given length n and a size k. This is technically 
possible because of the following 

Theorem: For any positive e the [n, k, d] CSS code as- 
sociated with randomly chosen classical codes C'2 C C\ C 



Z™ with dimC 2 = [2=*] and dimCi = [2±£] satisfies 
k/n>{l-e)R css {d/n) 

with a probability larger than 1 — 2~ n ( ae+ci ( n )) , where 
a is a positive constant independent of n and k. 

The situation is thus very much like the one in classical 
coding theory: randomly chosen subspaces of Z£ yield 
codes that asymptotically reach the Gilbert- Varshamov 
bound The theorem can be proven 0] along the 

same lines as the proof in 0] for the classical case, plus 
application of Mac Williams' theorem [l^j . 

Let an average over all CSS codes of length n and size k 
be defined via the uniform average over all pairs Ci C C\ 
of linear subspaces in with dimensions dimC2/i = 
[^S^]. According to the theorem, this average can be 
understood as an average over essentially all CSS codes of 
length n that asymptotically correct up to t qubit errors, 
where t = [^^-] is determined by k/n — R css (d/n). We 
denote this average by [•■■]"■ 

Applying the average to Eq. I|1U|) leads us finally to 



En 
w=t+ 



, OP™, (12) 
1 \w/ 



where we used the completeness relation a„„ = 1, 
and suppressed terms of order 2~'°^'. The physical in- 
terpretation of this expression is that more than t simul- 
taneous qubit errors cannot be corrected, and therefore 
coefficients cx vv with \v\ > t contribute to the residual 
error. Terms with non-diagonal coefficients a u ^ u ' turn 
out to be suppressed by a factor 2 t~ and therefore do 
not significantly contribute. Up to these exponentially 
small corrections the code averaged residual error is in- 
dependent of the encoded logical state I^l- 

In the following we consider [A]™ for a fixed ratio q = 
^ii. Inserting Eq. 10 into l|12|) immediately leads to 



[A]" 



E" ( n >y(i-*.r 

i — <w=qn \w/ 



(13) 

From this general formula for the residual error we can 
now draw conclusions on the performance of CSS codes. 

Independent errors correspond to the case of diverging 
distance r, where the decoherence parameter r r vanishes. 
In this limit the Gaussian in Eq. (|13fl shrinks to a nor- 
malized delta peak at x — 0. Hence, 



E 



r)p™(i- Po r- w , (w) 

-qn \<W J 



where p = (1 — e _r °)/2 is the error probability of a single 
spin. This expression, formerly derived by Calderbank 
and Shor [lfij . implies a binomially distributed residual 
error that decays exponentially with n as long as p < q. 
The point is that the constraint on the single-spin de- 
coherence Tq w 2p a < 2q is n- independent. In this 
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sense, CSS codes provide scalable quantum-error correc- 
tion with an exponentially small residual error |l5j . 

For finite r r the asymptotic value of the residual error 
can be easily derived from Eq. Ijl^fl by the observation 
that in the limit n — > oo the sum over w vanishes if p x < q 
and equals unity if p x > q. We obtain 



lim [A] n qn 



dx 



erfc 



Px>q 




where the last approximation is good for IV, To <C q. In 
sharp contrast to the independent errors discussed above 
for any finite T r (i.e., for any finite r) the residual error 
converges for large n to a finite constant (cf. Fig.^). 

According to quantum complexity theory 0] the total 
error in an m step quantum computation must be less 
than m -|0(i)l in order to produce useful results. If we 
therefore demand that the residual error must be limited 
by a A max (n) ~ b/n^ with some positive b and /i, we 
deduce from ifTBT) that T r must obey 



r r < T max (n) ~ qf (const + /xlnn) 



(16) 



This is no longer independent of n as for the uncorrelated 
model but approaches zero in the limit n — > oo, although 
slowly with the inverse logarithm of n. 

A further discussion must take into account that also 
T r by Eq. @ depends on n via the n-dependence of the 
maximum distance r and the observation time r. It is 
plausible to assume r oc n v with some exponent y > 1/3, 
and also r oc n v , because the observation time r must 
scale at least linearly with the time needed to transmit 
signals between qubits. We therefore set r = aro, r = 
aro, where a = (n/no) v , and write 

T,. = T(s, ar , ar , T, fi) = a ~ s r(s, r , To, aT, afl) . 

The second equality is obtained by rescaling the inte- 
gral J2J). Inspection of Eq. also reveals that for large 
effective temperatures aT r$ ,Tq the decoherence 
parameter is linear in temperature, and virtually inde- 
pendent on il. For large n condition (|16fl is therefore 
equivalent to 



r( s ,r ,r ,r)< const ^ ln „ (^)^- 2 ) 



(17) 



The inequality is violated if s < 2 and n is large. 
Consequently, in this regime the residual error exceeds 
A max (n). On the other hand, for s > 2 conditions i|17|) 
and (|l(jf) are satisfied for large n, meaning that here a 
sufficiently small residual error can be attained. We con- 
clude that for interaction paramter s > 2 scalable CSS 
error-correction is possible, while for s < 2 it is not. 

Apparently, noise correlations introduced by the spin- 
boson interaction are a strong hindrance for small values 
of s, whereas they are less harmful for larger s. An intu- 
itive explanation is that with increasing s spectral weight 
in J{uj) is shifted from lower to higher frequencies, and 



therefore the larger the value of s, the faster the disturb- 
ing noise correlations decay with distance. 

Note that the physical model we use is lacking dissipa- 
tive couplings, which corresponds to a systematic under- 
estimation of quantum noise. So, for a spin-boson model 
including dissipative couplings the residual error will be 
larger than the one calculated here. A second caveat 
is that we considered perfect quantum error-correction, 
while real systems will have to rely on fault-tolerant 
error-correcting schemes. It is not clear to what extend 
the results presented here do apply to those schemes. 
However, a naive guess might be based on the obser- 
vation that fault-tolerant quantum computation is, in 
some sense, simulating perfect quantum error-correction 
with imperfect quantum gates. This suggests that fault- 
tolerant schemes cannot better perform than perfect 
error-correction investigated here. Following these argu- 
ments, it is important to note that our results are not at 
odds with a recently proven threshold theorem [3| ■ Here 
the spin-boson interaction is unbounded and therefore 
does not belong to the class of interactions considered in 
Ref. 0. 

We wish to thank H. Moraal for numerous valuable 
discussions and M. R. Zirnbauer for critical reading of the 
manuscript. The work is supported by the SFB TR/12. 
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